We study the stability/establishment of dipolar mag- 
netostatic equilibrium configurations in new-born neutron 
stars (NSs) in dependence on the rotational velocity f2 
and on the initial angle a between rotation and magnetic 
axis. The NS is modeled as a sphere of a highly magne- 
tized (B r~j 10 15 G) incompressible fluid of uniform density 
which rotates rigidly. For the initial dipolar background 
magnetic field, which defines the magnetic axis, two dif- 
ferent configurations are assumed. We solve the 3D non- 
linear MHD equations by use of a spectral code. The prob- 
lem in dimensionless form is completely defined by the 
initial field strength (for a fixed field geometry), the mag- 
netic Prandtl number Pm, and the normalized rotation 
rate. The evolution of the magnetic and velocity fields 
is considered for initial magnetic field strengths charac- 
terized by the ratio of ohmic diffusion and initial Alfven 
■ travel times Tohm/'?"A.o ~ 1000, for Pm — 0.1,1,10, and 
, the ratio of rotation period and initial Alfven travel time, 
^ P/r A ,o = 0.012, 0.12, 1.2, 12 We find hints for the existence 
1—5 , of a unique stable dipolar magnetostatic configuration for 
\Q ' any specific a, independent of the initial field geometry. 

, Comparing NSs possessing the same field structure at the 
,— 1 ' end of their proto-NS phase, it turns out that sufficiently 
, fast rotating NSs (P ^ 6 ms) with a ^ 45° retain their 
magnetar field, while the others lose almost all of their ini- 
tial magnetic energy by transferring it into magnetic and 
\Q • kinetic energy of relatively small-scaled fields and con- 
, tinue their life as radio pulsars with a dipolar surface field 
)Q • of 10 12 ~ 13 G. 
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Abstract. 



1. Introduction 

The strength of a neutron star's (NS's) dipolar magnetic 
field, its inclination angle with respect to the NS's rota- 
tion axis, a, and the rotational velocity 17 are the quanti- 
ties determining the observability of pulsars as well as the 
power of their magneto-dipole radiation and/or the pul- 
sar wind, which e ventually set the spin-down behavior of 
an isolated pulsar l|Shapiro fe Teukolsk\ll98^1 . Therefore, 
the determination of these (more or less exactly observ- 
able) quantities in their temporal evolution is subject of 
scientific debates since the discovery of pulsars. 

From the observational point of view there exists only 
little information about the magnetic inclination a im- 
mediately after birth of a NS. This is owed to the poor 
observability of new-born N Ss, which are surrou nded 
by opaque ma terial. Recently, lHarding et al.1 |2p03) and 
iHardina <)2006(l investigated the emission geometry of very 
young and energetic pulsars and found that five out of six 
have a < 30°. The 1600 years old pulsar J1119-6127, e.g., 
has an inclination angle a = 19° and a dipolar magnetic 
field strength of 4.2-10 13 G. Age and dipolar magnetic field 
strength correspond to the observed spin period of 0.42 s, 
provided the field has been constant during the lifetime of 
the NS. Given that constancy an initial spin p eriod of less 
than a millisecond can be estimated (see, e.g.. lCok>i et alJ 

Unfortunately, there is even less observational evidence 
about a in magnetars because their slow rotation impedes 
the detection of open field line regions. The typically wide 
X-ray pulses come likely from a very wide beam. If the 
X-ray emission region is not too far above the surface and 
the pulse peak corresponds to the magnetic pole, then 
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that wid e beam may indicate that the inclination angle a 
is small l|Zhandl2006)) . 

In order to understand those, yet vague, observational 
evidences of relatively small values of a for magnetars and 
young energetic pulsars, and to get at all an idea about 
the "initial" magnetic field structure of NSs, the study of 
the establishment/stability of magnctostatic equilibrium 
configurations immediately after a NS's birth is necessary. 
In this very early phase, when the core is not yet super- 
fluid, consists of npe-matter, and when the crust is com- 
pletely liquid, an equilibrium configuration may indeed be 
reached. It defines the initial magnetic field structure and 
strength, especially that of the dipolar component, which 
plays a decisive role for spin-down behavior and observ- 
ability of the NS. 

In new-born NSs, the rotation velocity exceeds the 
hydromagnetic (or Alfven) one, ua, in large regions, even 
for magnetic field strengths in the order of magnetar fields 
(B < 10 15 G). Therefore, the effect which rotation may 
have on the establishment of a magnetostatic equilibrium 
in young NSs has to be studied in detail. 

It was IWrightl lll973h . who postulated the stabil- 
ity of a magnetic field configuration the poloidal and 
toroidal constituents of which are of comparable strength. 
This conclusion was s uppor ted by subsequent studies, 
e.g..lMarkev fc Tavlerl l|l973j) or iBraithwaite fc Nord lund 
(200(J. Tayler and co-authors performed a series of ana- 
lytical stability analyses of magnetic fields in stars : lTavlerl 
lll973l) considered the stability of a purely toroidal ax- 
isymmetric field without rotation and found that it is 
always unstable with respect to azimuthally large-scale 
(m = 1) perturbations. The growth rate of this instabil- 
ity is on the order of the Alfven frequency — VA/r, 
where va = B/^Airp, B the magnetic field strength, p 
the density and r the radial coordinate. This "Tayler- 
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instability " , meanwhile basic ingredient of the "Spruit- 
dynamo" acts locally in meridional 
(r, 9) planes but g lobally in t he azimu t hal ((f ) direction. 
iMarkev fc Tavlerl (|l9T3^ and I Wright! l)l973f) performed 
similar adiabatic stability analyses for an axisymmetric 
poloidal field and concluded its general instability if at 
least some of its field lines are closed within the star. 
The instability generates small scale field structures on 
the sam e time scale as for a toroi dal field which dissipate 
rapidly. IMarkev fc Tavl^l ^197.^ discuss the possibility 
that sufficiently rapid rotation and/or the simultaneous 
existence of equipollent toroidal and poloidal components 
ma y stabilize the axisymm etric background fields. 

Pitts fc Tavlerl <)l985|) explored in ideal MHD and 
cylindrical geometry the effect of rigid rotation on the sta- 
bility of a field consisting of a toroidal component which 
increases with the distance from the magnetic axis, and/or 
a uniform z (i.e., poloidal) component for different ana- 
lytically feasible special cases: For a = both the purely 
toroidal field and the combined toroidal-poloidal one show 
qualitatively the same stabilizing effect of the rotation. For 
a = 90°, however, a purely toroidal field reveals a much 
less stabilizing effect of the rotation. They studied fur- 
ther the stabilizing effects of a stable stratification of the 
star's density profile in the case a = 0° and for a purely 
toroidal field which may incr ease non-line arly with axial 
distance. As already shown bv lTavlerl l|l957!) . such a field is 
stable even in the absence of rotation for sufficient small 
wavelengths of a perturbation along the magnetic axis. 
Larger scaled perturbations, however, cause instabilities 
which can be count eracted by sufficiently rapid rotation. 
On the other hand, IPitts fc Tavlerl l|l985^ found that at 
the same rotation rate, which stabilizes the original insta- 
bilities, rotation introduces a new instability. Therefore, 
they concluded, that "rotation is unlikely to lead to com- 
plete stability of general magnetic field configurations" . 

In 2005 the interest in the stability of magnetostatic 
equili bria in stars became r evitali ze d by a serie s of pa - 
pers: iBraithwaite fc Spruit! ll2006(k iBraithwaitel l|2005|) : 
iBraithwaite fc Nordlundl l|2006|) . The latter authors stud- 
ied the stability of magnetic fields in the radiative in- 
teriors of Ap stars and found, that any random field 
is generally unstable but evolves in a stably stratified 
star towards a "twisted torus" configuration with approx- 
imately equipollent toroida l and poloidal components. 
IBraithwaite fc Spruit! l)200fil) considered the stability of 
MHD configurations in ultra-magnetized NSs, the mag- 
netars. Although they found that, quite similar to the Ap 
stars, stable magnetic fields exist when being concentrated 
to a relatively small region around the center of the star 
and after having evolved into a poloidal-toroidal torus 
shape, these fields cannot have magnetar strength unless 
their initial strength exceeds 10 18 G. This estimate re- 
lies on two field-diminishing eff ects in the simulations of 
IBraithwaite fc Nordlundl l)2006|) : Firstly, a decrease of typ- 
ically < 10 appears during the relaxation into the torus- 
like structure; secondly, these fields are strongly concen- 
trated in the central region and decay by a factor of up to 



1000 towards the surface. However, NSs born with mag- 
netic fields of ~ 10 18 G are unlikely. On the other hand 
an internal field of (only) magnetar strength would need 
a much too long time to diffuse to the surface. 

Until now the effect of rotation on the stability of 
purely poloidal field configuratio ns has not be e n con - 
sidered in gen eral. The results of IPitts fc Tavlerl ljl985j) . 
iLuvtenl l|l988j) and lCurrv fc Pudritd l)l995|) do not apply 
directly since these works consider combined toroidal and 
poloidal fields , the latter being homog eneous. In the sole 
case in which Iciirrv fc Pudritd l)l995|) consider an inho- 
mogeneous poloidal field, the Tayler instability was not 
found. The work of Braith waite and Spr uit doesn't in- 
clude rotation except in IBraithwaite! l)2005h who considers 
toroidal fields only. 

Since rotation in new-born NSs can for sure not be 
neglected and the dipolar magnetic field determines both 
observability and rotational evolution, we will study the 
consequences of rotation for the stability of such fields. In 
the next section we describe and justify our model assump- 
tions, present the set of equations to be solved together 
with the initial conditions for the models under considera- 
tion and sketch our numerical approach. After presenting 
the results for the different MHD configurations in Sect.O 
we finally discuss the results. 

2. Model, basic equations, numerical method 

2.1. Description of the model 

We study the stability of magnetostatic configurations in 
a model applicable to strongly magnetized, rotating new- 
born NSs which have just established their (almost) sta- 
tionary density profile according to the thermodynamic 
properties of NS matter (expressed by an equation of 
state). It is a highly conductive normal-fluid npe liquid 
considered incompressible, a perhaps not too crude ap- 
proximation for young NSs. Furthermore, we assume a 
constant density profile and rigid rotation throughout the 
star, that is, all internal motions to be damped out ini- 
tially. 

Let us discuss these model assumptions in more detail: 

1 . We assume that the crust is still completely liquid and 
that the core matter has not yet performed the phase 
transition into the superfluid state. The onset of super- 
fluidity and/or crystallization limits the period during 
which a magnetostatic equilibrium can establish undis- 
turbedly to at most a few hours. Although the temper- 
ature of the transition into the superfluid phase, T c , is 
rather uncertain, it probably lies in wide regions of the 
NS between 10 9 and 10 10 K. According to t he approxi- 
mate a nalytical formulae given in Sect. 3 of lPage et al! 
(|2005() . the period after which the NS has cooled down 
to these temperatures is 17 s ... 0.5 yrs. On the other 
hand, a characteristic Alfven time is on the order of 
ta ~ 0.02 ... 20 s, depending on the background field 
strength. Wisely assuming that T c ^ 8 • 10 9 K for the 
majority of non-exotic equations of state, we believe 



4 U. Geppert and M. Rheinhardt: Magnetars versus Radio Pulsars 



that there is ample time for the MHD configuration to 
reach either a stable state or to go through an unstable 
phase dissipating magnetic energy very efficiently (see 
iBraithwaite & Spruitj l20061 . The same is valid for the 
onset of crystallization at the crust-core interface. A 
typical value for the melting temperature at the den- 
sity there, p w 1.5 • 10 14 g cm" 3 , is T m w 8 ■ 10 9 K 
ilonesl 11990) so that about the same cooling period 
may elapse until the onset of crystallization in the crust 
as passes until the transition into the superfluid state 
in the core and inner crust. Hence, using again the 
above mentioned formulae, for at least 1000 seconds 
after the NS's birth the state of its matter is that of a 
normal-fluid npe liquid. 

2. At the end of the PNS phase a very strong magnetic 
field is present i n the star, either due to PNS dynamo 
action (see, e.g.. lRheinhardt fc Geppertll2005t and ref- 
erences therein) or simply by flux conservation during 
the collapse. The field may comprise components of 
very different scale lengths, but is likely to possess a 
strong, large scale, say, dipolar fraction which defines 
the magnetic axis of the star. 

3. It is very likely, too, that the NS is rapidly rotating 
immediately after its birth. Since braking effects had 
not yet time to act efficiently, it follows from angular 
momentum conservation (see, e.g., model calcula tions 
of the progenitor's rotation in iHeeer et al.ll2005)) that 
the initial spin period of th e NS lies in the range of 
1 ... 10 ms. lOtt et alJ l|2005|) studied the evolution of 
rotational profiles in proto-NSs taking into account 
magnetic effects in parameterized form. They found 
that a proto-NS rotating differentially at the begin- 
ning will reach very fast the state of rigid rotation 
by angular momentum redistribution due to turbu- 
lent (magnetically induced) viscous dissipation, so that 
at the end of the proto-NS phase, when the definite, 
non-convective, configuration of the new-bo rn NS is 
reach ed, it rotates most likely rigidly (see also lShapirol 
12000ft . 

4. In order to let the assumption of incompressibility be 
justified the fluid velocity should be well below the 
sound velocity. We estimate the latter to be in the 
range c s — 8 • 10 8 cm s _1 . . . < speed of light for the 
density rang i ng fro m 10 12 to 3T0 16 g cm~ 3 according to 
iBavm et al 1 lll97ll) . Of course, the results of our model 
have to be checked a posteriori to be consistent with 
this requirement. The same has to be done with respect 
to the criterion that the propagation time of sound 
waves across the star must be short in comparison with 
the time scales of the velocity. 

5. The assumption of constant density is surely less well 
justi fied. However, for non-exotic equations of state 
(see IBavm et al.lll97l|) the change in density is only 
about one order of magnitude within three quarters 
of the NS volume. Again a-posteriori, it has to be 
checked whether the bulk of the flow is concentrated 
within that core volume. As a future refinement of the 



model, realistic density profiles should be included in 
the framework of the anelastic approximation. 
6. As both the magnetic and kinetic energies are small in 
comparison with the gravitational binding, any devi- 
ation of the NS from sphericity is neglected. Another 
approximation consists in the assumption of a perfect 
vacuum in the NS's surroundings which implies neces- 
sarily a vanishing radial component of the velocity at 
the surface. Both approximations together imply that 
only two of the three components of the force on a 
surface element can be required to be zero (see Eq. 

2.2. Basic equations 

We decided to solve the non-linear MHD equations as pre- 
sented below although a linear analysis would be sufficient 
to detect instabilities. Hence, we solve an initial-value 
problem where the background field together with some 
arbitrary perturbations form the initial condition. Since 
we imprint only perturbations whose energy is very small 
in comparison with that of the background field, there is 
always a linear stage from which the same growth/decay 
rates can be derived as from a stability analysis based 
on the linearized equations. The advantage of solving the 
non-linear equations consists in the additional informa- 
tions we obtain with respect to the duration of the lin- 
ear stage and the saturation of the "final" MHD config- 
uration which may represent another, now stable, equi- 
librium. As the symmetry of the fastest growing eigen- 
modes of the perturba tions is a-priori unknown (cf., e.g., 
iMarkev fc Tavlerlll973h we have to solve in 3D without 
any symmetry constraints. 

The preceding considerations result in the following 
ruling MHD equations for magnetic field B and velocity 
u and corresponding boundary conditions written in di- 
mcnsionless form: 



dB 



dt 

curlB = 
div.B = 
[B] = 



— AB + curl(w x B) 



for 
for 



r < 1 



r > 1 

Vr 

for r = 1 



(1) 



du 



1 



-V(p- ^Q 2 r 2 sin 2 tf + $) 

(uV)it + Pm Ait 
2£le z x u + curl.B x B 



> for r < 1 (2) 



div u = 



u r = 
(D(u) ■ e r )# tV = 



for r = 1 



(3) 



Here, length is normalized on the NS radius R, time on the 
magnetic diffusion time Tohm = R 2 /v> where r\ — c 2 /4ira 
is the magnetic diffusivity and a the electric conductivity. 
[.] denotes the jump of a quantity across a surface. The 
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magnetic field is measured in units of (^■np) 1 l 2 r\j R, the 
velocity in units oi-q/R, the pressure p in units of pr] 2 /R 2 , 
and the gravitational potential $ in units of rj 2 / R 2 . D(u) 

denotes the deformation tensor, (r, #, (p) are spherical co- 
ordinates, the polar axis — 0, ft e 2 ) of which coincides 
with the axis of rotation. 

Equations Q and J5J) contain two parameters which 
together with the initial field strength Bq (for a fixed field 
geometry) define the problem completely: the magnetic 
Prandtl number Pm and the normalized rotation rate O. 

Pm, is the ratio of kinematic viscosity, v, and magnetic 
diffusivity, 77, and is chosen to be 0.1, 1, and 10, resp., what 
is partly dictated by numerical restrictions, but on the 
other hand represents a subset of the NS-relevant range 
determined by the prevailing densities (p = (1...2.8) • 
10 14 gem -3 ) and temperatures (~ 10 11 - 9 K). These values 
result in Pm = 0.1. . 6 ■ 10 s llCutler fc Lindblomlll987t 
lYakovlev fc Shalvbkovlll99l reflecting mainly the strong 
temperature dependence of Pm oc T~ 4 . 

As we rely on a dimensionless representation, the two 
remaining parameters are best expressed as ratios of char- 
acteristic times: f2 by qp = P/tj^ q, P the rotation pe- 
riod of the NS, and the initial field amplitude B by 
QB = TA,o/min{r^ cay ,r| ecay }, where r A ,o = R/v A ,o 
is the Alfven crossing time related to the initial field. 
Note, that the "decay times", r^ e , r£ e , are differ- 
ent from the often used magnetic and viscous "diffusion 
times", Tohm = R? ' jf], TVisc = R 2 /v, by factors l/V 2 and 
0.06676, respectively. In the following we use Tdecay for 

mm {^fccay: T dccay}- 

For qs Q we take values ~ 10~ 2 , which are of course by 
far too large for highly magnetized young NSs where this 
ratio can go down to 10~ 16 for a surface magnetic field 
of 10 15 G. However, as long as the growth/decay times 
of the examined perturbations of the initial state are in 
the order of at most a few Alfven times, dissipation of the 
background state does not affect their linear stage (quasi- 
stationary approximation) . 

The artificially enhanced dissipation is a con- 
cession to numerical feasibility only, in order to 
avoid excessive requiremen ts for spatial resolution (cf. 
iBraithwaite fc Spruit! 1200^3 . where qB Q was chosen to be 
0.1). Hence, the use of the results from the nonlinear stage 
has to be considered with caution, (see Sect.0J). 

For all models in the present study we choose the fol- 
lowing parameters 

R = 10 6 cm , p = 2 • 10 14 g/cm 3 , B(t = 0) = 10 15 G, (4) 

where B is the de-normalized r.m.s. value. Therefore, 
r A,o — 0-05 s. With rotation periods chosen between 0.6 
ms and 0.6 s, where the former value is close to the gener- 
ally accepted minimum for new-born NSs, but the latter 
is in that respect not very likely, these model parameters 
result in qp = 0.012 ... 12. As a reference and for com- 
parison with known results for non-rotating stars we also 
consider qp — 00. 



2.3. Initial conditions 

Each initial condition consists of an axisymmetric dipolar 
background field the stability of which is examined, and 
imprinted magnetic perturbations. The fluid is assumed 
to be at rest initially. The initial magnetic inclination a is 
defined by the dipole axis. Of course, for a ^ the back- 
ground fields shown in Figs.nand[21have to be considered 
properly rotated. The energy of the perturbations is set to 
0.1% of the magnetic background energy, and their cho- 
sen geometry ensures that the initial state has no preferred 
equatorial symmetry. Their azimuthal dependence is char- 
acterized by m — 6, 7. This apparently arbitrary choice is 
motivated by numerical feasibility but at least doesn't af- 
fect the validity of the results for the linear stage. 

2.3.1. Internal uniform field 




Fig. 1. Internal uniform field according to JSJ in a merid- 
ional plane. The arrow indicated below marks a magnetic 
field of 10 15 G. 

As a first, more academic example we consider a field 
uniform throughout the sphere, but being a dipolar vac- 
uum field outside: 

B = B e z for r < 1 

/ cos \ 

B = -Bo curl r x V ( - 2 j for r > 1 (5) 

[B r ] =0 for r = 1 

Here, Bq denotes the polar surface field strength. Of 
course, the continuity of the normal component of B has 
to be required, but the tangential components remain dis- 
continuous and give rise to surface currents. We have to 
concede, however, that the used numerical method en- 
forces continuity of all components of B on the surface 
whereby we in fact used a field which in a thin shell imme- 
diately below the surface deviates from homogeneity (see 
Fig.^). We include this model field into our study because 
it was the one consideration s on magnetic stability i n NSs 
were first ex emplified on bylFlowers fc Rudermanl l|l977h 
and because [Braithwaite fc Spruitl l|2006f) reported on its 
instability in the non-rotating case. 
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2.3.2. Dipolar magnetostatic equilibrium field 




t = 0.0000 



Fig. 2. Dipolar magnetostatic equilibrium field according 
to (JSJ) . The arrow indicated below marks a magnetic field 
of 1.71 • 10 15 G. 



A more realistic initial configuration consists in the 
poloidal magnetostatic equilibrium field the angular de- 
pendence of which is the same as f or a dipolar field 
iRohertJTflSlllRneinfwdt et aliEool see Fig.©: 



B = B a curlr x V ( ^(3r 3 — 5r) cost? ) for r < 1 



B = -Bn curlr x V 



cos $ 
~2^~ 



for r > 1 



IB] 







for r = 1 
(6) 

As the Lorentz force of this field is a pure gradient it can 
in the constant -density case be balanced by the quan- 
tity -V(p - {fl 2 r 2 sin 2 i9)/2 + $) (see Eq. ©). If stable, 
this state could be a final equilibrium to which an arbi- 
trary initial configuration relaxes during the fluid stage of 
a new ly born NS. As known f rom analytical studies, how- 
ever, llMarkev & Tavlerlll97.^ any purely poloidal field is 
in the absence of rotation surely unstable, hence the evo- 
lution of the field 10 under the influence of rotation is a 
crucial question. 

2.4. Numerical method 

Because the star is modeled as a spherical body and 
incompressibility is assumed, the problem is especially 
suited to be tackled by a spectral method. The one used 
here employs an expansion of B and u into their modes 
of free decay Bi(r), Ui(r). The Bi(r) are the time- 
independent parts of the eigensolutions of with u = 0. 
Similarly, the Ui (r) are obtained as the time-independent 
parts of the eigensolutions of J2J, © with B = 0, ft = 0, 
and the term («V)ti canceled. It can be shown, that both 
sets of modes form orthogonal systems of vector func- 
tions, which are complete with respect to all vector func- 
tions obeying the respect ive boundary cond itions and be- 
ing solenoidal (see, e.g., iRheinhardtl 1 1 99W . That is, the 



solutions can be represented as the infinite sums 

oo 

B(r,t) = J2 b i{t)Bi(r) , 

i=l 
oo 

u(r,t) = y^Uj(t)uj(r) 



(7) 



with 



AB, + XfBi = , divB 4 = 
AUi + XfUi = Vpi , div u t = 



(8) 



where for numerical feasibility the sums have to be 
suitably truncated. The expansion functions are either 
poloidal ("P") or toroidal ("T"), each of them rep- 
resentable by means of a single spherical harmonic 



B 



Tm 



B 



Pm 



N^cuTlrxV Jn (,4fr)Y, 



= -Ki rxVj n (n u J t)y; 



iiT , 



(9) 



U 



71=1, 



^ p curlrxV(j n (/i«f r)- 
; m = 0, . . . , n ; 1 = 1,. 



Here, j n (x) is the nth order spherical Bessel function of 
the first kind; the eigenvalues /^ p , /-i^ T , y^J , , ordered 
increasingly with the index I, are determined such that the 
boundary conditions are obeyed. The collective index i in 
0) is in © resolved into {n, m, I, pi}, where pt stands for 
"P" or "T" . The quantities are normalization factors. 

The resulting system of ordinary differential equations 
for the mode amplitudes bi(t) and Ui(t) is integrated by 
use of a Runge-Kutta scheme of 4th order with adaptive 
step-size control. Usually, the truncation of the expansion 
series ensured that the first neglected modes contained 
not more than 0.1 % of the energy of the dominating one 
and that the energy spectra showed a clear cutoff, mostly 
being exponential. 

3. Results 

3.1. Solution characteristics 

Equations CQl~© allow of solutions with special symme- 
tries with respect to the equatorial plane and the rotation 
axis, resp.: By the symbol "A" we denote those solutions 
for which B/u is equatorially antisymmetric/symmetric. 
Analogously, "S" denotes solutions with both fields being 
equatorially symmetric. Referring to the azimuthal depen- 
dence, the fields can either be both axisymmetric or super- 
positions of constituents each characterized by a <p depen- 
dence ~ exp(imt/?). In the general case all m > occur, 
but there are special solutions for which only constituents 
for a subset of m values are non-zero. In the extreme case 
that this subset consists only of one m, the fields are de- 
noted by this value or, if there is an additional equatorial 
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symmetry, by "Am" or "Sm" . For axisymmctric (with re- 
spect to the rotation axis) background fields the extreme 
case emerges just for the linearized versions of Eqs. Q 
and (J2J- That is, each of the corresponding eigensolutions 
can be represented by a single m. 

In order to characterize the solutions more generally we 
use apart from energies and spectra the following quanti- 
ties related to symmetry properties of the fields. The "par- 
ity" II reflects symmetry with respect to the equatorial 
plane: II = ±1 means perfect symmetry /antisymmetry. 
The deviation from axisymmetry is described by the "non- 
axisymmetry" M : M — 0, 1 refers to perfect axisymmetry 
and missing of any axisymmetric constituent, respectively. 
The general definitions are: 

E — E 4.- E° 

JT — s y m -^anti M — 1 — 

E sym + E ant i E 

Here, E synl , E ant {, E°, and E are the energies of the sym- 
metric, antisymmetric and axisymmetric (m = 0) parts of 
a field and its total energy, respectively. 

The runs were usually extended in time to the mag- 
netic decay time, that is, to t£, . For the calculation of 
ta we use the r.m.s. value of the magnetic field. In order 
to clarify the difference between the purely ohmic decay 
of a field and its evolution under the influence of the flows 
driven by its own Lorentz force, we relate the final mag- 
netic energies of both cases after equal evolution times, 
where the final energy of the purely ohmic case is denoted 
by -Emag 1 - Growth/decay times for a certain perturbation 
symmetry (to) extracted from the linear stages (see Tables 
12 and |3J) are those of the fastest growing of all modes for 
this to. 



the final magnetic energy is almost equal to the one re- 
sulting from ohmic decay only, are stable as confirmed by 
examination of the linear stage. All other cases turn out 
to be unstable ones. Note, that in case of extremely fast 
rotation (qp = 0.012 or P = 0.6 ms) even background 
fields with a — 90° can be stabilized. 

3.2.2. Dependence on the magnetic inclination 

If we fix gp=0.12,1.2 and Pm = 1, but consider a = 45° 
and a = 90° it turns out that at the same rotation rate 
and increasing a the stabilizing effect of rotation vanishes; 
for gp=0.12, a = 45° still instability is observed, but ap- 
parently close to its marginal value. 

3.2.3. Dependence on the magnetic Prandtl number 

Since both diffusivities are heavily overestimated in our 
calculations, we decided to adjust the normalization such 
that in comparison to the case Pm = 1 the case Pm = 0.1 
corresponds to a reduction of the viscosity by a factor 10 
while for Pm = 10 the same holds for the magnetic diffu- 
sivity. Thus we reduced in both cases the total dissipation 
in comparison with the case Pm = 1. Consequently, af- 
ter the same elapsed period of time the total energy for 
Pm =/= 1 is always larger than for Pm = 1. Reducing Pm 
from 1 to 0.1 obviously doesn't cause qualitative changes 
of the stability properties. An increase of Pm, however, 
impedes stabilization, as visible from the cases with a = 0, 
qp = 0.12. As Pm actually increases with progressing cool- 
ing, perhaps an even more rapid rotation is necessary to 
ensure stabilization in the cooler stages. 



3.2. Dipolar magnetostatic equilibrium 

In this model the strength of the initial background field, 
B , is chosen such 1 that 

fl.5-10- 2 for Pto=1,10 



As the symmetry parameters II and M are defined with 
respect to the rotation axis, they depend on a. For 
a = 0°,45°,90° the initial values II = -1,0,1 and 
M = 0, 0.5, 1 are obtained, respectively. The results for 
Pm = 0.1,1,10 and the various a are summarized in 
Table [T] 

3.2.1. Dependence on the rotation period 

For fixed Pm and a the final magnetic energy increases 
and the final kinetic energy decreases with decreasing P. 
Simultaneously the magnetic symmetry parameters ap- 
proach those of the initial field. The three cases, in which 

1 For Pm = 1 and as r A ,o/ min{r^ cay , r| ecay } = ^/Wp ■ 
r\ max(7r 2 , Pm -15) the normalization field is then 10 -3 • 
B T .m. B .(t — 0) so that the initial dimensionless field is 10 3 . 



3.2.4. A stabilized case 

Let us consider the case qp = 0.12, Pm = 1, a — which 
for the model parameters Q corresponds to P = 6 ms, 
that is, a very fast spinning new-born NS. Of course, the 
non-linear evolution is in general influenced by the arbi- 
trariness of the initial conditions. However, in the present 
case it is without relevance since all perturbations, apart 
from a short initial "switch-on" , decay very fast as shown 
in Fig. |3J There, time-changing energy spectra £™ (i), 
(t) are shown in which all modes with the same to (cf. 
Sect. 12.4)) are summed up to yield a single spectral contri- 
bution labeled by m. Note, that the temporary growth of 
the kinetic energy for m — can be attributed to the error 
introduced by representing the background equilibrium as 
a truncated series of modes. 

Due to the axisymmetry of the background field and 
all coefficients in JJJ and J2J|, each eigenmode of the lin- 
earized problem is with respect to its dependence on tp 
uniquely characterized by a single m. Therefore, our spec- 
tral representation with respect to m allows identifying 
the decay time of the slowest decaying eigenmode for each 
m from the data of the linear stage. When performing 
this, one has to skip an initial phase in which the decay 
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Table 1. Results for the dipolar magnetostatic equilibrium model. All quantities are taken after a period of T^ e with 
some exceptions denoted in the footnotes. Subscripts "kin" and "mag" refer to velocity and magnetic-field related 
quantities, respectively. Boldface symmetry symbols indicate, that the geometries of the initial and final magnetic 
configurations coincide thus giving a further hint for stability. The calculations for the non-rotating NS, qp = oo, 
were performed only for a = because in this case the choice of the axis is of course arbitrary. The results shown for 
shorter integration times, T cn d/T^ c < 1, allow nevertheless unambiguous conclusions about stability. 



Pm a(°) 


QP 


' ni a,g / ' ni ag 


mag 


n 

mag 


Illcin 


M 

mag 


Mkin 


B symmetry 




OO 


0.0011 


6.22 


0.44 


0.19 


0.63 


0.86 


mixed 




12. 


0.0024 


1.93 


-0.44 


-0.15 


0.17 


0.85 


mixed 





1.2 


0.072 


0.013 


_i 


0.66 


0.049 


0.84 


~ AO 


0.1 


0.12 


0.98 


0.0003 


-1 


1 





0.9 


AO 


-t'j 


1.2 






n Qu 




u.uo 




144-1-A.tJU 




12. 


0.0014 


6.51 


0.009 


0.1 


0.43 


0.47 


mixed 


90 


1.2 


0.0342 


0.0033 


0.85 


0.4 


0.8 


0.97 


mixed 




0.12 


0.139 


0.0009 


1 


-0.38 


1 


0.94 


SI 




oo 


0.0002 


4.67 


-0.06 


0.6 


0.4 


0.99 


mixed 




12. 


0.0004 


6.6 


-0.5 


0.8 


0.2 


0.88 


mixed 





1.2 


0.0076 


0.0007 


-0.99 


0.66 


0.06 


0.91 


mixed 




0.12 


n qs 

u.yo 


n nnnfn 

u.uuuuo 


_l 


I 


n 

u 


n ooq 

u.uuy 


AO 


1 45 


1.2 


0.075 


0.001 


-0.46 


-0.14 


0.27 


0.97 


mixed 




0.12 


0.824 


0.00005 


-0.24 


0.64 


0.38 


0.64 


mixed 




12. 


0.0033 


0.64 


0.5 


0.78 


0.77 


1 


mixed 




1.2 


0.043 


0.003 


0.88 


0.67 


0.81 


0.65 


mixed 


90 


0.12 


0.14 


0.00014 


1 


1 


1 


0.59 


SI 




0.012 


0.98 


0.0013 


1 


0.97 


1 


1 


SI 




oo 


0.00097 


0.305 


0.155 


-0.512 


0.86 


0.95 


mixed 





1.2 


0.028 


0.0041 


-0.40 


-0.32 


0.41 


0.91 


mixed 


10 


0.12 a 


0.77 


0.0017 


-0.95 


0.045 


0.053 


0.94 


~ AO 




1.2 


0.033 


0.0033 


0.19 


0.23 


0.68 


0.93 


mixed 


90 


0.12 6 


0.74 


0.0017 


0.99 


0.82 


0.997 


0.906 


~ AO 



Tend — 0.72r^ C a 
Tend — 0.37~ dcca y 



is not purely exponential. This is a consequence of the 
arbitrary initial conditions for the perturbations feeding 
several eigenmodes with different decay times simultane- 
ously. However, there is always an instant of time after 
which the slowest decaying mode survives as the only one 
yielding a single, well-defined decay time. 

The decay times are shown in Tabled in units of ta q. 
The value for m = 1 is not included as apparently the 
slowest decaying mode belonging to m = 1 is not yet 
clearly dominating over the faster decaying m — 1 modes. 

Moreover, Fig. [21 shows the field parities and non- 
axisymmetries as functions of time. Clearly, the magnetic 
field retains its symmetry completely whereas the velocity 
starts with equatorial antisymmetry and complete non- 
axisymmetry as can be explained from the interaction 
of the perturbations with the background field, but ap- 
proaches very fast its final equatorially symmetric, ax- 
isymmetric geometry. Note, that even the maximum of 
iiW is much smaller than E mag . A comparison of the ini- 
tial and the final state reveals that the final state coincides 



almost exactly with the one resulting from free decay of 
the initial state alone. Hence, in all we can safely conclude 
that the considered dipolar equilibrium is indeed stabilized 
by rotation. 

Table 2. Decay times Ty of the slowest decaying eigen- 
modes depending on m for the dipolar equilibrium model 
with qp = 0.12, a = and Pm = 1 (a stabilized case). 



m2345678 
40.97 24.73 20.82 13.94 10.89 8.50 8.31 



3.2.5. An unstable case 

It is instructive to compare the former case with a much 
slower rotating NS: P/ta,o = 12, i.e. the same setup as 
before, but P = 0.6 s. 
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Fig. 3. Temporal evolution of the disturbed dipolar magnetostatic equilibrium for qp=0.12, Pm = 1, a = 0. Time is 
in units of rohm, energy in units of prj 2 /R 2 . Subscripts "kin" and "mag" refer to velocity and magnetic-field related 
quantities, respectively. Upper two rows: magnetic and kinetic energies, each further subdivided in its poloidal and 
toroidal part; third row: m spectra of the magnetic and velocity fields, E™ ag , E^ n , respectively. Lower left and right 
panels: parity and non-axisymmetry. 



Inspecting Fig. it becomes clear that the slower ro- 
tation is no longer capable of stabilizing the background 
field. The kinetic energy of the evolving flow is in its max- 
imum about 10 times larger than in the stable case with 



rapid rotation. In the temporal evolution of the energies 
we can distinguish three different phases: initially, up to 
« 6ta,0) the effect of the perturbations is negligible and 
the slow ohmic decay dominates. Afterwards, the rapidly 
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increasing perturbations drain out energy from the back- 
ground field very efficiently. From s» 8ta,o on, saturation 
sets in followed by a decelerating decay of all energy con- 
stituents to low values. 

Checking the validity of the incompressibility assump- 
tion, we find that, except for a very thin outer shell where 
in real NSs the density would be less than 10 8 g cm -3 , 
the maximum velocities are well below the speed of sound. 
However, we have to concede that for realistic viscosities 
this assumption is likely to be violated during the fully de- 
veloped non-linear stage, at least locally. A compressible 
treatment would then be necessary. On the other hand, 
during the early phase of the non-linear stage, that is, 
before the kinetic energy peaks, the observed time scale 
is short enough to be unaffected by the assumed unrealis- 
tically high diffusivities. Therefore, we are confident that 
up to this moment our results are reliable with respect to 
magnitudes, too. 

The time-changing m-spectra in Fig.0]show nearly ex- 
ponential growth of the perturbations on the linear stage 
from which growth times were deduced (see Table yj. 
While for low m the behavior of / ta,o with respect to m 
is non-monotonous, a continuous decrease is observed be- 
yond m = 6. Thi s is in accordance w ith predictions from 
analytical work (jPitts fc Tavlerlll985|) . 

The time dependences of parity and non-axisymmetry 
shown in Fig. ^demonstrate the increasing deviation from 
the equilibrium, too. 

Table 3. Growth times of the fastest growing eigen- 
modes depending on m for the dipolar equilibrium model 
with qp = 12, a = and Pm = 1 (an unstable case). 



m 


1 


2 


3 


4 5 


6 


7 


8 


1~A,0 


0.45 


0.25 


0.17 


0.23 0.48 


0.90 


0.79 


0.32 



effectively than by ohmic diffusion only, thereby deceler- 
ating the overall decay. 

A comparison of both the final spectra and field ge- 
ometry of the result for a = 0, qp = 0.12 with the cor- 
responding result for the dipolar equilibrium model gives 
strong evidence for a tendency to approach the same state 
for t — > oo (see Fig. |HJ|. Not only the difference in the ini- 
tial geometry (cf. Figs. and EJ, but also the difference 
in the initial energies is obviously equalized after having 
gone through the non-linear stage. This is confirmed by 
the similarity of the spectra of the two results the remain- 
ing differences of which disappear in the further course of 
time. Fig. shows the final field for the internal uniform 
field model which is almost identical with the final field of 
the poloidal equilibrium model. The relative r.m.s. value 
of the difference of both fields is only 1.6 %. The same 
coincidence is found for a = 45°, 90°, qp — 0.012. 




| = 7.10e+02 



Fig. 6. Final field geometry for the internal uniform field 
model in a meridional plane (the field is almost exactly 
axisymmetric) for a = 0, qp = 0.12, Pm = 1. With 
the denormalization based on Eq. the maximum field 
strength is about 7.1 • 10 14 G. 



3.3. Internal uniform field 

Here, Bq is chosen such that 

7"A,o/ T dccay = 8 ■ 10~ 3 

and its symmetry parameters n and M coincide with those 
of the dipolar equilibrium field. The results for Pm = 1 
and a = 0°, 45°, 90° are summarized in Table 0] The gen- 
eral tendencies are the same as for the dipolar equilibrium 
field. Again, the significant influence of a on the stability 
is proven. In the cases a = 0, qp = 0.12 and a = 45°, 
qp = 0.012 the final magnetic energy is slightly larger in 
the MHD case in comparison with the purely ohmic decay. 
This somewhat surprising result is not in contradiction to 
energy conservation because it cannot be excluded that 
the flow driven by the Lorentz force reorganizes the mag- 
netic field in such a way that small scale components just 
below the surface (cf. Sect.EQJ ar e smoothed out more 



4. Discussion and Conclusions 

The results of the present study suggest the following ma- 
jor conclusion: Never again in an isolated NS's life its mag- 
netic field may suffer so drastic changes as during the first 
at most 1000 seconds after birth, provided the dipolar field 
present immediately after the proto-NS phase has mag- 
netar strength. Depending on the rotational velocity and 
the inclination of the field with respect to the rotation 
axis, the inborn field will either survive that period as a 
magnetar field or lose almost all of its initial magnetic en- 
ergy by transferring it into magnetic and kinetic energy 
of quite small-scaled fields eventually being consumed by 
ohmic and viscous dissipation. After the end of the un- 
stable phase, large-scaled magnetic fields dominate again, 
but have a much lower strength. 

We found that purely dipolar fields in highly magne- 
tized new-born NSs are unstable on a time-scale of a few 
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Table 4. Results for the internal uniform field model. Only Pm = 1 was considered. For explanations see Tabled 
Boldface symmetry symbols denote here that the final geometry coincides with the one of the dipolar equilibrium. 



a ( ) 


Qp 


P 1 pOhm 
^mag / J-^mag 


^kin / u mag 


TT 

11 mag 


TT, 


M 

J" mag 


M, ■ 

J" km 


-D by LLlLLLChLy 





oo 
12. 
1.2 
0.12 


0.00009 
0.001 
0.083 
1.02 


147.6 
2.92 
0.0034 
0.00002 


-0.41 
-0.82 

-1 

-1 


-0.43 
-0.44 
0.22 
1 


0.82 
0.11 
0.03 



0.81 
0.86 
0.93 
0.0001 


mixed 
mixed 
mixed 
AO 


45 


0.12 
0.012 


0.56 
1.006 


0.00006 
0.0023 


-0.33 
-0.011 


-0.15 
-1 


0.34 
0.49 


0.92 
1 


mixed 
A0+S1 


90 


0.12 
0.012 


0.193 
0.976 


0.0005 
0.0013 


1 

1 


1 
1 


1 
1 


0.91 
0.997 


SI 
SI 



Alfven times unless fast rotation, i.e., flj^/Cl 0.1, acts 
stabilizing. The latter condition, however, is only a neces- 
sary one, fulfilled for NSs which have at the end of their 
proto-NS phase a poloidal field strength of ~ 10 15 G and 
a rotation period of 1 ... 6 ms. 

As a sufficient condition, the above presented results 
yield that rotation is only efficiently stabilizing as long 
as the inclination angle between rotation and magnetic 
axis a < 45°. Therefore, we may speculate that those NSs 
whose inborn field is not too much bent away from the 
rotation axis and which rotate with a initial period ^ 6ms 
will continue their life as magnetars, while for quite similar 
NSs, which, however, are born as orthogonal rotators, even 
rotation periods « 6ms will not stabilize the poloidal field 
against its substantial reduction down to standard pulsar 
strengths. As another conclusion we suggest that for each 
inclination angle a there exists a "final" , apart from ohmic 
decay, stable field configuration, whatever the inborn field 
shape has been. Indeed, they seem to be identical up to a 
rot ation by a. 

iFlowers fc Rudermanl (jl97?f ) consider qualitatively the 
relaxation of a uniform field into its lowest energetic state. 
Assuming a relatively large Alfven time scale of ~ 100 s 
(in 1977 magnetar field strengths were inconceivable!) and 
a rotation period of the new-born NS « 6 ms, they con- 
clude that the time for the readjustment of that field is 
about 10 days, clearly beyond the onset of crystallization 
which impedes the relaxation already after a few 1000 sec- 
onds. Their relaxation time scale is (except for a factor of 
2) that of the Tayler instability, ta, enhanced by a factor 
ta IP due to the strong effect of the Coriolis force I Spruit! 
1999). Making use of this argumentation with respect to 
magnetar field strengths, however, this estimate for the 
relaxation time scale reduces to ~ 1 s, i.e., the initially 
uniform magnetic field would have ample time to estab- 
lish its energetically lowest state which is supposed to be 
no longer dominated by a large scale field. In contrast, our 
study shows that this will not happen for a rotation period 
of < 6 ms unless a ^ 45°. For an understanding of this dis- 
crepancy the beha vior of the so-calle d Levitron™ can be 
in structive (cf. e.g.lJones et alJll99^ : Similar to the setup 
of IFlowers fc Rudermanl (jl97?f ) it consists of two perma- 
nent, mutually repelling magnets where a stable (levita- 



tion) situation can be achieved only via rotation. However, 
it seems questionable at all whether the analogy with per- 
manent magnets is appropriate for NSs because the latter 
can be described in terms of ideal MHD characterized by 
a frozen-in field whereas the former are characterized by 
frozen-in currents. 

In this paper we consider the stability of purely 
poloidal fields. However, it is very unlikely, that in na- 
ture such a configuration exists. More likely is that 
the poloidal field is accompanied by a toroidal one, 
which, if it is of comparable strength, may act stabi- 
lizing as well llTavlerl Il980t IBraithwaite fc SpTuil 120041: 
iBraithwaite k, Nordlundl l2006[) . The dependence of the 
stability of such a configuration on rotation and inclina- 
tion will be studied in a forthcoming paper. 

It is well-known, that the effect of a NS's stably strat- 
ified d£nsity_p_rofiIe exerts another stabilizing effect (see, 
e.g., IPitts fc Tavlerlll985h . This is due to the strong sup- 
pression of motions in the direction of gravity being an 
additional constraint for the development of an instabil- 
ity. Hence, the neglect of this effect in the present study 
justifies the more our conclusions about the stabilizing 
action of rapid rotation. This could, however, imply that 
the marginal P, which discriminates between stability and 
instability, shifts to somewhat larger values. In principle, 
compressibility, which is the cause for stratification, affects 
any MHD process and may give rise to specific instabil- 
ities, e.g. the buoyancy (Parker) one, not existent in an 
incompressible setup. Thus, deviation from the assump- 
tion of uniform density may add both a stabilizing and 
a destabilizing tendency. Although we are, guided by our 
estimate of the ratio of sound speed to fluid velocity, quite 
sure that the stabilizing effect prevails, a quantitative cer- 
tainty can only be gained on the basis of a numerical study 
with an (at least) anelastic model. 

Apart from the magnetic field the vigorous flows oc- 
curring in the unstable situations may result in observ- 
able consequences: When extrapolating the maximum ve- 
locities found in the simulations to lower, more realistic 
viscosities, it seems likely that in some regions of the NS 
velocities of *S 10 8 cms -1 can be reached. Supposed that 
the initial perturbation for the instability considered in 
this paper is given by the (decaying) PNS convection, it 
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is conceivable that the geometry of the latter is somehow 
"transferred" to the one of the Lorentz-force induced flow. 
Moreover, as the rotation rates of the PNS and the newly 
born NS are almost identical, it is allowed to speculate 
that similar to the PNS convection the demonstrated un- 
stable decay of a strong magnetic field may be a source 
of gravitational waves. Applying likewise the estimates of 
Miill er fc Jankal lll997l) . we conclude that for a duration 
from seconds to a few minutes a detectable gravitational 
wave signal can be expected rising at about 20. . . 30 s after 
bounce. 

Our study is based on the assumption, that newly born 
NSs inherit from their proto-NS phase surface magnetic 
fields > 10 15 G. This state can be reached either by flux 
conservation if the progenitor was already highly mag- 
netized, by a PNS dynamo, or by other, still unknown, 
processes. 

If the NSs are born with a strong dipolar field then 

1. those whose rotation period is less than ~ 6ms and 
whose magnetic inclination angle a is smaller than 
~ 45° will retain their extremely large surface field 
strength and appear, after a rapid spin down, as mag- 
netars; 

2. those which rotate less rapidly, say with P ^ 6ms 
and/or for which is a <; 45° will lose almost all of their 
inborn magnetic energy and appear as radio pulsars. 

It turns out that rotation is likely to be the only stabilizing 
agent which allows of the existence of magnetars whereas 
the st able configurations found by iBraithwaite fe Spruit! 
l)2006(l are less suited to support such strong surface fields. 

In the light of the present study, the much smaller 
number of observed magnetars in comparison with that 
of all the other observed realizations of NSs seems to be 
a consequence of the fact that only a small fraction of 
all new-born NSs are rotating as fast as or faster than 
P ~ 6ms. Of course, this conclusion is still premature, 
lacks broad observational support, does not consider other 
selection effects, and has to be confirmed by the study of 
both more realistic NS models and a denser mesh of initial 
a and P values as well as by studies of the evolutionary 
stages 'supernova progenitor' - 'supernova' - 'proto-NS'. 
These considerations should definitely include the effects 
of magnetic fields and return a better idea about the dis- 
tribution of magnetic field strengths, rotation periods and 
inclination angles at the beginning of a fully evolved NS's 
life. 

Let us now relax the assumption that neutron stars 
have at the end of the PNS phase a field of ^ 10 15 G 
and adopt instead fields at that moment which are two 
to three orders of magnitude smaller. When additionally 
conceding that the initial rotation period could be as large 
as, say, 60 ms, the ratio qp would be still at most ~ 10~ 2 . 
As our results demonstrate, such fields would be always 
stable irrespective of the magnetic inclination a. 
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Fig. 4. Temporal evolution of the disturbed dipolar magnetostatic equilibrium. qp=12 1 Pm = 1, a = 0. For further 
explanations see Fig. |2| Note that the numerical resolution with respect to m in the vicinity of the kinetic energy 
peak, though being worse than the usually obeyed standard (cf. end of Sect. 12.4(1 . still ensures an exponential spectral 
decay of ~ two orders of magnitude. Using the parameters Q t = 0.01 corresponds to 0.5 seconds! 
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Fig. 5. Field geometries of magnetic field (left) and flow (right) around the end of the dramatic field reduction phase 
(see upper panels of Fig.0]at t ps 0.01) in a meridional plane for the dipolar equilibrium model with gp=12, Pm = 1, 
a = 0. Arrows indicate vector components parallel to the paper plane. Their maxima are 6.93 • 10 14 G and 8.3 • 10 6 
cm s -1 , respectively. Colors encode the field moduli: the brightest tone corresponds to 7.2 • 10 14 G and 1.2 • 10 7 cm 
s _1 , respectively. 



